Phonon spectrum, thermal expansion and heat capacity of UO2 from first-principles 
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We report first-principles calculations of the phonon dispersion spectrum, thermal expansion, 
and heat capacity of uranium dioxide. The so-called direct method, based on the quasiharmonic 
approximation, is used to calculate the phonon frequencies within a density functional framework for 
the electronic structure. The phonon dispersions calculated at the theoretical equilibrium volume 
agree well with experimental dispersions. The computed phonon density of states (DOS) compare 
reasonably well with measurement data, as do also the calculated frequencies of the Raman and 
infrared active modes including the LO/TO splitting. To study the pressure dependence of the 
phonon frequencies we calculate phonon dispersions for several lattice constants. Our computed 
phonon spectra demonstrate the opening of a gap between the optical and acoustic modes induced 
by pressure. Taking into account the phonon contribution to the total free energy of UO2 its thermal 
expansion coefficient and heat capacity have been ab initio computed. Both quantities are in good 
agreement with available experimental data for temperatures up to about 500 K. 

PACS numbers: 63.20.dk, 65.40.Ba, 65.40.De 



I. INTRODUCTION 

Over the last few decades UO2 has been one of most 
widely studied actinide oxides due to its technological 
importance as standard fuel material used in nuclear re- 
actors. There exists currently considerable interest in un- 
derstanding the behavior of nuclear fuel in reactors which 
is a complex phenomenon, influenced by a large num- 
ber of materials' properties, such as thermomechanical 
strength, chemical stability, microstructure, and defects. 
Especially, knowledge of the fuel's thermodynamic prop- 
erties, such as specific heat, thermal expansion, and ther- 
mal conductivity, is essential to evaluate the fuel's per- 
formance in nuclear reactors-li^ These thermodynamic 
quantities are directly related to the lattice dynamics of 
the fuel materiaLl^ 

Dolling et alWweie the first to measure phonon disper- 
sion curves of UO2 , using the inelastic neutron scattering 
technique in 1965; their seminal article has become the 
standard reference for uranium dioxide's phonon spec- 
trum. Later the vibrational properties of UO2 were in- 
vestigated in detail by Schoenes, lousing infrared and Ra- 
man spectroscopic techniques. A good agreement with 
phonon frequencies obtained from inelastic neutron scat- 
tering was observed!^ More recently, Livneh and Sterei^I' 
studied the influence of pressure on the Raman scattering 
in UO2 and LivnelP^ demonstrated the resonant coupling 
between longitudinal optical (LO) phonons and U'*"'' crys- 
tal field excitations in a Raman spectroscopic investiga- 
tion. A theoretical investigation of the phonon spectra of 
UO2 was reported recently by Yin and Savraso\ff^ who 
employed a combination of a density-functional-theory 
(DFT) based technique and a many-body approach. Ac- 
cording to their results, the low thermal conductivity 
of UO2 stems from the large anharmonicity of the LO 



modes resulting in no cont ributi on from these modes in 
the heat transfer. Goel et a/.'^^'^ investigated the phonon 
properties of UO2 using an empirical interatomic po- 
tential based on the shell model and observed that the 
calculated thermodynamic properties including the spe- 
cific heat are in good agreement with available exper- 
imental data. Devey^^ employed recently the general- 
ized gradient approximation with additional Coulomb U 
(GGA-f J7) to compute the main phonon mode frequen- 
cies at the Brillouin zone center which were in reason- 
able agreement with experimental data. Very recently, 
Sanati et al}'^ used the GGA and GGA-I-C/ approaches 
to investigate phonon density of states and elastic and 
thermal constants, which were found to be in reason- 
ably good agreement with experimental data. In spite of 
the already performed studies, further investigations are 
needed. Especially, the full dispersions of the phonons 
in reciprocal space have not yet been considered. Also, 
important quantities such as the thermal expansion coef- 
ficient and heat capacity are directly related to the lattice 
vibrations but these quantities have not yet been studied 
ah initio from the calculated phonon spectrum. 

The objective of this study is to contribute to a detailed 
understanding of the lattice vibrations of UO2. Using 
the first-principles approach, based on the DFT we have 
calculated phonon dispersion curves and phonon density 
of states of UO2. The calculated phonon properties are 
compared with the available experimental data from in- 
elastic neutron scattering and Raman spectroscopy along 
with a detailed discussion. Furthermore, several thermo- 
dynamic properties have been computed taking the influ- 
ence of lattice vibrations into account. Here, we report 
the lattice contribution to the heat capacity as function 
of temperature as well as temperature and volume (in the 
quasiharmonic approximation). The dependence of the 
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total free energy on the lattice constant of UO2 as a func- 
tion of temperature has calculated, from which we derive 
the thermal expansion coefficient. The thermal expan- 
sion coefficient as well as lattice heat capacity compare 
favorably to available experimental data up to 500 K, 
which is the temperature range in which the influence of 
anharmonicity can be neglected. 



II. COMPUTATIONAL METHODOLOGY 

The electronic structure of UO2 has been discussed in 
the past years .'32H2?DFT calculations within the general- 
ized gradient approximation (GGA) underestimate the 
influence of the strong on-site Coulomb repulsion be- 
tween the 5/ electrons. An improved 5/ electronic struc- 
ture description can be obtained with the GGA+U ap- 
proach, in which a supplementary on-site Coulomb repul- 
sion term is added; this approach correctly gives the elec- 
tronic band gap of UOz-E^ISlIllHlol while the GGA+U ap- 
proach would appear preferable for description of UO2 's 
electronic structure, we encountered specific problems 
when using this method. Some of the phonon branches 
became negative away from the zone center. This artifact 
might be related to the occupation matrix of 5/ states 
that would require an additional stabilizing constraint 
in the GGA+U method.'^^H! Using conversely the spin- 
polarized GGA approach, we found that such difhcul- 
ties did not occur. The phonon dispersion spectrum pre- 
sented below is hence computed with the GGA exchange- 
correlation for antiferromagnetically ordered UO2 and is 
found to be in good agreement with experiment .^^ 

Here, we have determined the phonon dispersion 
curves and density of states (DOS) in t he qu asiharmonic 
approximation using the direct method.'^^El gy displac- 
ing one atom in a supercell (of 96 atoms) from its equilib- 
rium position, non-vanishing Hellmann-Feynman forces 
were generated. Due to the high symmetry of the face- 
centered cubic (fee) lattice of UO2, only one atom for 
uranium (U) and for oxygen (O) was needed to be dis- 
placed. The actual shift of the atoms in the supercell had 
an amplitude of 0.03 A and was taken along the [001] di- 
rection only, on account of the cubic symmetry of UO2. 
In the calculation of the resulting forces we employed 
the projector augmented wave (PAW) pseudopotential 
approach within the Vienna Ab-initi o Sim ulation Pack- 
age (VASP).l2MIlThe PHONON codS^ has been used 
to extract the force constant matrix from the Hellmann- 
Feynman forces and to subsequently calculate the phonon 
dispersion curves and DOS. 

For the thermodynamic quantities we consider the to- 
tal free energy of UO2, including the phonon contribu- 
tion, 

F{e, T) = U{e) + i^P''°"(e, T) + F^'(e, T), (1) 

where -F(e, T) is the Helmholtz free energy at a given 
strain e. The phonon free energy contribution i^p''°" is 



expressed as 

/•OO 

FP'"'"{e,T) = / dug{uj,e)[nuj/2 

+fcBTln(l-e-'""/^'^^)], (2) 

where g{u!,e) is the phonon DOS, computed as mentioned 
above. We note that the free electronic energy, F'^\e, T), 
is not considered in the present study, because the ther- 
mal electronic contribution is known to be negligible in 
the temperature ran ge up to 1000 K, which is the range of 
interest in this work.l^^l^Sl 'pjjg static lattice energy [/(e) 
appearing in Eq. ([I]) can be expressed as 

C/(e) = C/o + F^C„-e,e„ (3) 

where Uq is the static lattice energy at zero strain, Cij 
are the elastic constants, and V is the equilibrium volume 
at T = K. The static lattice energ ies have also been 
calculated using the VASP code\^^ 

In our calculations we have used a 2x2x2 su- 
percell containing 96 atoms with a 4x4x4 fc-point 
mesh in the Brillouin zone (BZ). The Perdew-Wang 
parametrization"^^ of the GGA functional was used. The 
kinetic energy cut-off for the plane waves was set at 600 
eV and the energy criterion used for convergence was 
10^^ eV. The force acting on each ion was converged until 
less than 0.01 eV/A. Once the phonon DOS has been cal- 
culated, the thermal expansion of UO2 can be evaluated 
straightforwardly. First, the phonon DOS with static lat- 
tice energy is calculated for several volumes around the 
T — K equilibrium volume. Subsequently, the total 
free energies are calculated for these different volumes at 
constant temperature using Eqs. ([l])-(|3|. After the free 
energy has been calculated its minimum gives the corre- 
sponding equilibrium volume at the considered tempera- 
ture. By repeating the process for different temperatures, 
the thermal expansion coefficient a defined by 

"(^) = adT = 3ydT' 

is obtained; here a is the lattice constant. 

A further thermodynamic quantity, the lattice contri- 
bution to the specific heat can be derived from ^'^^^'^^ at 
a fixed temperature in the quasiharmonic approximation. 

III. RESULTS AND DISCUSSION 

A. Ground-state properties of UO2 

UO2 crystallizes in the cubic fluorite structure (CaF2) 
belonging to the Fm3m space group (no. 225) and there 
are three atoms per primitive unit cell with one U atom 
(Wyckoff position 4a) and two inequivalent O atoms 
(Wyckoff positions 8c). Therefore, there are generally 



TABLE I: Calculated equilibrium lattice constant a (in A) 
and bulk modulus B (in GPa) of UO2. Theoretical values 
obtained in this work are compared to values from molecu- 
lar dynamics simulations'^ (MD), as well as to experimental 
datEpaEH (Exp.). 





a (A) 


B (GPa) 


This work 
MD (300 K) 
Exp. (300 K) 


5.406 
5.47*1 

5 4';J39|40| 


184 
IS*' 

192F'l98,E2l20'PSl 



nine phonon branches. Before turning to the descrip- 
tion and analysis of the calculated lattice dynamics let us 
briefly consider the ground-state properties of UO2. As 
mentioned above the employed DFT framework is that 
of the spin-polarized GGA. The calculated equilibrium 
lattice constant a and bulk modulus 5, which we have 
obtained by a Birch-Murnaghan 3'"'' order fitj^ are pre- 
sented in Table |lj where these are compared to experi- 
mental lattice properties .'S^'^ Our calculated equilibrium 
volume as well as the bulk modulus compare reasonably 
well with the experimental dat a a nd with results from 
molecular dynamics simulations.'^ Compared to the ex- 
perimental lattice constant the lattice constant computed 
here is 1.2% smaller. This can be attributed to a too 
strong binding of the 5/ orbitals which become too much 
delocalized in the spin-polarized-GGA approach. In the 
GGA+C/ approach the 5/ orbitals are more localized and 
their contribution to the bonding reduced/^ ^"^ which 
leads to a theoretical lattice parameter which is larger 
than the experimental one.'^S 



B. Phonon spectrum of UO2 

Figure [T] shows the measured and ah initio calculated 
phonon dispersion curves of UO2 along high-symmetry 
points of the fee BZ. The black lines are the calculated 
phonon dispersions obtained for the different lattice pa- 
rameters, the blue open symbols are the experimentally 
measured data from inelastic neutron scattering^ and 
red closed symbols are data from Raman scattering.li^l 
Note that M denotes the additional point (1, 1, 0) that 
was included by Dolling et al.^ but that is outside of 
the first BZ. The calculated phonon dispersions along 
the high-symmetry lines X-W-F have been added to il- 
lustrate the whole dispersion of phonon states in the BZ. 
The three branches in the low frequency region are the 
transverse acoustic (TA) and longitudinal acoustic (LA) 
modes that belong to vibrations of the U atom with its 
relatively heavy mass. The other branches are optical 
modes with higher frequencies that are mainly associated 
with lattice vibrations of O atoms and can be labeled as 
TOl, LOl and T02 and L02, respectively (see Fig.[T]), as 
there are two inequivalent O atoms in the unit cell. The 
three panels in Figure [T] illustrate the volume dependence 
of the phonon frequencies, which have been calculated at 
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FIG. 1: (Color online) Calculated (full lines) and measured 
phonon dispersion curves (open symbols, from Ref. |9| along 
high-symmetry directions in the fee Brillouin zone. The mea- 
sured Raman and infrared modes (Ref. I12p at the F point are 
depicted by solid (red) circles. The notation of the special 
points is M: (1,1,0), F: (0,0,0), X: (1,0,0), W: (1, |, 0), 
and L: (|, \, |). 



three different lattice constants, the experimental one, 
a — 5.47 A, the optimized theoretical one, a — 5.406 A, 
and a = 5.33 A, a selected smaller lattice constant. 

In the surrounding of the F point the agreement be- 
tween the calculated and measured phonon frequencies 
is very good. A first discrepancy between measured 
and computed dispersions is observed along the M— F 
high-symmetry direction. Our calculations predict three 
acoustic branches whereas in experiment there are two 
branches. It could be that a splitting of the low-lying 
TA branch near the M point could not be sufficiently 
resolved in the experiment. This might be related to 
the inelastic neutron scattering technique, which might 
be affected by the size of samples and as well as a rela- 
tively low neutron flux (see, e.g., Refs. 43 and 44). Also, 
as mentioned by Dolling et al. the frequency measure- 
ments of the phonons with the neutron technique might 
be impeded in the zone boundary regions. The LO and 
TO branches agree reasonably well with experiment, in 
particular for the optimized theoretical lattice parame- 
ter. There are some discrepancies in the positions of the 
branches at the zone boundary X and at the M point. 
One of the TOl branches turning up from the F point 
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to the M point has not been detected in the experiment. 
Along F— L the agreement is quite good. The top-most 
(L02) branch deviates most between calculation and the 
experiment. This high-lying branch has the largest ex- 
perimental uncertainty. Nonetheless, the ab initio cal- 
culated branch has more dispersion than present in the 
measurement and, except for the zone center, it falls out- 
side of the experimental error bar. 

The phonon frequencies calculated at a = 5.406 A and 
a = 5.47 A are very similar to each other and agree rea- 
sonably well with the measurement data. A notable dif- 
ference between the two sets of dispersion curves appears 
however in the frequency gap at the zone boundaries of 
M and X. When the lattice constant decreases with pres- 
sure from a = 5.47 A to a = 5.33 A a pressure-induced 
phonon softening occurs. The frequency gap between 
LA and TOl modes at the zone boundary X point and 
at the M point is increased as the pressure increases. At 
a = 5.47 A the LA and TOl modes almost approach 
each other, whereas at a = 5.406 A, a gap is predicted to 
exist between the LA and TOl modes. A small or van- 
ishing gap between these modes is in accordance with the 
measurement. We also note that with increased lattice 
constant the negative slope of the LOl branch is remark- 
ably increased along the F— X symmetry line. This find- 
ing suggests that the propagation of the LOl phonon is 
significantly restrained as the lattice constant increases. 

In Fig. [l] we furthermore depicted by the red solid 
circles at the F point the frequ encies obtained by Ra- 
man and infrared measurement s.fiSEll Using group the- 
ory analysis these active Raman and infrared modes can 
be decomposed into irreducible representations of the 
(O^) point group, as lT2g + 2Ti„. The U atom con- 
tributes only to the infrared active mode (Ti„), whereas 
the O atom contributes to both, the infrared and Raman 
mode, Ti„ and T2g, respectively. Both these modes are 
triply degenerate. The frequencies of these modes are 
summarized in Table [iTj together wit h result s from MD 
simulations^^ and experimental results.l^ESIil! Xo account 
for the LO/TO splitting the Lyddane-Sachs-Tellei'^ re- 
lation was used on the basis of effective charges. These 
were chosen to be 3 for the uranium atom and —1.5 for 
oxygen. These values are consistent with values used in 
the literature.'^ The macroscopic electric field in UO2 
splits the infrared-active optical modes into TO and LO 
components. Frequencies of TO modes are calculated in 
a straightforward manner within the direct method but 
the LO modes can only be obtained via introduction of 
a non-analytical terrrPSl into the dynamical matrix. In 
general, this term depends on the Born effective charge 
tensor and the electronic part of the dielectric function 
(high-frequency dielectric constant). Details regarding 
this procedure can be found elsewhere, see e.g. Ref. SHI 

The agreement of the computed Raman optical mode 
as well as the infrared TO mode with available experi- 
mental data is quite good, see Table |ll| The calculated 
frequency of the infrared LO mode is somewhat smaller 
than the experimental results, indicating that the effec- 



TABLE 11: The frequencies of the Raman (T29) and infra-red 
(Tin) vibrational modes of UO2 at the F-point. Theoretical 
results (this work) are given for three lattice constants a and 
compared to frequencies obtained from molecular dynamics 
simulatiorP (MD) as well as experimental (exp.) frequencies 
(given in THz). 



Mode 




this work 


MD^^ 


exp. 


a (A) 


5.33 


5.406 


5.47 


5.47 


5.47 


Ti„(TO) 


9.27 


8.49 


7.78 


7.62 


8.3^, 8.m 










8 5^^ 


Ti„(LO) 


16.5 


15.82 


15.26 


17.28 


16.6#a, 16.713 










17.3^, nM 


T29 


14.22 


13.4 


12.80 


14.04 


12.9^33, I3.4P 



five charge and dielectric constant taken from literature 
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do not fully account for the correct — = In ad- 

dition there is a spread of ca. 0.7 THz in the m easured 
values of the LO infrared mode frequencies! ^"*"'''' ! This is, 
however, not sufficient to account for the underestima- 
tion of the experimental values by about 2 THz that is 
observed for a — 5.470 A (upper panel of Fig. [l]). Using 
the calculated theoretical volume (middle panel of Fig. 
[ij this deviation is reduced by 25%. A similar error of ca. 
1.5 THz for the ah initio calculated infrared frequencies 
of another insulating compound, CsNiFa, was reported 
in a recent study. 

Figure [2] shows the calculated and measured phonon 
DOS. The blue, red, and green lines with square, tri- 
angle, and diamond symbols indicate the phonon DOS 
computed for a=5.330 A, 5.406 A, and 5.470 A, respec- 
tively. The experimental datcP (a « 5.47 A, T=296K) 
are plotted with the black line. The U contribution to 
the calculated phonon DOSs gives rise to a higher in- 
tensity and narrower peak widths in the lower frequency 
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"a =5.406 A 

— a= 5.470 A 

— experiment 




Frequency (THz) 

FIG. 2: (Color online) The theoretical phonon density of 
states (DOS) of UO2 computed for three lattice constants 
a = 5.330, 5.406, and 5.406 A, compared to the measured^ 
DOS (for a « 5.47 A, at r=296 K). 
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region. The more broadened DOS with lower intensity 
that occurs in the higher frequency region is mainly de- 
rived from the oxygen atoms. A notable difference be- 
tween the phonon DOS at the three lattice parameters 
is the size and position of the phonon gap occurring for 
frequencies of about 6 — 7 THz. For the experimental 
lattice parameter a = 5.47 A the computed gap practi- 
cally closes. The experimental phonon DOS spectrum at 
this lattice parameter shows a minimum at about 6 THz, 
in reasonable agreement, considering some experimental 
broadening. We note that the trend of decreasing gap 
with larger lattice constant continues, leading to a clos- 
ing of the gap computed for larger lattice constants (not 
shown). Overall, the computed phonons DOS of both 
the theoretical equilibrium (a = 5.406 A) and the experi- 
mental lattice parameter are in good accordance with the 
measured spectrum. The phonon DOS of the theoretical 
lattice parameter agrees best with the experimental data 
at higher frequencies (7 to 13 THz), where the peaks co- 
incide with the measured ones. As mentioned earlier, the 
L02 mode lies both lower in the computed spectra and 
is more dispersive than in the measurements. We note 
that the recent GGA-I-J7 calculations^^ provide a sharper 
DOS peak at 17 THz, due to a flatter L02 dispersion near 
the zone boundaries. At the zone center the L02 branch 
lies however much deeper than in the experiment, at 10 
THz (vs. 17 THz in experiment). 
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FIG. 3: (Color online) Top: calculated and experimental ther- 
mal expansion coefficient a{T) of UO2. The experimental 
data are those of Taylor!^ Bottom: Computed temperature- 
dependent lattice parameter of UO2 (open squares) and spline 
fit function. 



[5T]) , which may add an additional influence on the lattice 
parameter. 



D. Thermodynamic properties of UO2 



C. Thermal expansion of UO2 

The calculated phonon DOS enables us to evaluate 
some thermodynamic quantities which depend on the lat- 
tice vibrations. We start with the thermal expansion. 
The phonon contribution to the total free energy of UO2 
increases with increasing temperature and hence becomes 
progressively responsible for changes of the lattice pa- 
rameters. To compute the thermal expansion of UO2 
we have first computed the total free energy, including 
the phonon contribution, for various lattice parameters, 
from which we computed the temperature-dependent lat- 
tice constant. Figure [s] (bottom) shows the calculated 
variation of the lattice constant with temperature. The 
red curve gives the spline interpolation of the calculated 
lattice constants, shown by the symbols. The thermal ex- 
pansion coefficient a{T) — a^^da/dT was subsequently 
evaluated by differentiating the spline fit. The upper 
panel of Fig. [3] shows the calculated thermal expansion 
coefficient, which is in good agreement with experimental 
datal^" up to 500 K. The deviation between the calculated 
and measured data slightly increases above 500 K and be- 
comes significant at around 1000 K. This might be due 
to an increased electronic contribution to the thermal ex- 
pansion. At very low temperatures, in the region between 
and 50 K, a deviation is also observed between the cal- 
culated and measured data. The origin of this deviation 
is not unambiguously clear. We note however that UO2 
undergoes a magnetic phase transition at 31 K (see Ref. 



Next, we employ the Helmholtz free energy to compute 
the lattice heat capacity at constant volume (Cy) and at 
constant pressure (Cp). Fig. |4] shows theoretical results 
for Cy, computed with the harmonic approximation, and 
Cp, computed with the qua siharm onic approximation, 
as well as experimental result J^^Ull for Cp up to 1000 K. 
First the lattice contribution to the specific heat Cv was 
computed as a function of temperature for the equilib- 
rium lattice constant, a = 5.406 A. Subsequently, the 
specific heat at constant pressure was derived according 
to 

Cp = Cv + 9a^Ba^T, (5) 

where a{T), a, and B are the calculated linear thermal 
expansion coefficient, the equilibrium lattice constant, 
and the bulk modulus (see, e.g., Ref. [55]). Fig. |4] il- 
lustrates that there is a very good agreement between 
the cornputed heat capacity Cp and the experimental 
data.'S^HSll Note that the sharp anomaly in the experi- 
mental data at T « 31 K is due to the aforementioned 
magnetic phase transition,!^!' which effect is not included 
in the calculations. Clearly, the specific heat at constant 
pressure is in much better agreement with the experi- 
mental data at higher temperatures than Cy, which is 
mainly due to the thermal expansion of UO2. Conversely, 
evaluating Cp for the theoretical equilibrium lattice con- 
stant or for the experimental lattice constant a = 5.47 
A only gives very minor differences. The computed Cp 
curves fall somewhat below the experimental Cp data 



6 



90 



70 r 
"5 60 

E 

3-50 



40 
CO 
Q. 

1 20 



- 




■ 

■ ■ ,* 






- • 
















■ Cp , exp. Gronvold 
• Cp, exp. Huntzicker 




r 






C^ (a = 5.406 A) 

--. Cp 






/jjr 




. _ . , Cp using exp. a 




- • 


If 








- • 1 










- • ,g 
' • JT 

•/ 












1 , , 1 , 


, 1 , , 1 , , 1 , 


1 < < 1 < . 1 . 





100 200 300 400 500 600 700 800 900 1000 

Temperature (K) 

FIG. 4: (Color online) The lattice contribution to the heat ca- 
pacity of UO2 at constant volume, Cv , computed within the 
harmonic approximation (full curve), and a constant pres- 
sure, Cp, computed with the quasiharmonic approximation 
(dashed lines; for a = 5.406 A and the experimental a = 5.47 
A). The experimental data for lower and higher temperatures, 
full circles and full squares, are taken from Huntzicker and 
Westrum (Ref.l52p and Gr0nvold et al. (Ref. 53), respectively. 



for temperatures in the range of 400 to 1000 K. The ura- 
nium atoms in UO2 are in the paramagnetic state at 
higher temperatures .1121411 Therefore, the remaining dif- 
ference between experimental and computed data above 
T > 400 K could be due to the magnetic entropy con- 
tribution to the specific heat or, alternatively, it could 
be due to the anharmonic effects. The magnetic en- 
tropy contribution to the specific heat was investigated 
for T = 200 to 300 K in Refs. 51 and 53, At higher 
temperatures it saturates to approximately S = i? In3, 
corroborating a magnetic J = 1 state on the uranium 
atoms. It provides to a relatively small magnetic entropy 
contribution that would lead to a small increase of the 
computed Cp data (by about 3 Jmol^^K^^). 



IV. CONCLUSIONS 

We have performed first-principles calculations to in- 
vestigate the lattice vibrations and their contribution 



to thermal properties of UO2. We find that the cal- 
culated phonon dispersions are in good agreement with 
experimental dispersion^ measured using inelastic neu- 
tron scattering. Computing the phonon dispersions for 
various lattice constants, we observed a softening of the 
phonon frequencies with decreasing lattice constant. Fur- 
thermore, the band gap between TA and LO modes at 
high-symmetry zone-boundary points are found to de- 
pend significantly on the volume. This gap almost closes 
at a = 5.47 A, consistent with a pseudogap detected in 
the inelastic neutron experiment. Also, Raman and in- 
frared active modes have been determined as a function of 
volume. The agreement with experimental data and with 
results obtained with molecular dynamics simulations is 
overall very good, with an exception of the infrared LO 
mode that appears underestimated in our first-principles 
calculations. Including the phonon contribution to the 
free energy, the heat capacity and the thermal expansion 
coefficient of UO2 have been computed. Both thermal 
quantities are found to agree well with experimental data 
for temperatures up to 500 K. The good correspondence 
of the computed and measured thermal data exemplifies 
the feasibility of performing first-principles modeling of 
the thermal properties of the important nuclear fuel ma- 
terial UO2. 
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